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Abstract. The cosmic rays differential intensity inside the heliosphere, for energy below 30 
GeV/nuc, depends on solar activity and interplanetary magnetic field polarity. This variation, 
termed solar modulation, is described using a 2-D (radius and colatitude) Monte Carlo approach 
for solving the Parker transport equation that includes diffusion, convection, magnetic drift 
and adiabatic energy loss. Since the whole transport is strongly related to the interplanetary 
magnetic field (IMF) structure, a better understanding of his description is needed in order 
to reproduce the cosmic rays intensity at the Earth, as well as outside the ecliptic plane. In 
this work an interplanetary magnetic field model including the standard description on ecliptic 
region and a polar correction is presented. This treatment of the IMF, implemented in the 
HelMod Monte Carlo code (version 2.0), was used to determine the effects on the differential 
intensity of Proton at 1 AU and allowed one to investigate how latitudinal gradients of proton 
intensities, observed in the inner heliosphere with the Ulysses spacecraft during 1995, can be 
affected by the modification of the IMF in the polar regions. 



1. Introduction 

The Solar Modulation, due to the solar activity, affects the Local Interstellar Spectrum (LIS) 
of Galactic Cosmic Rays (GCR) typically at energies lower than 30 GeV/nucl. This process, 
described by means of the Parker equation (e.g., see [H [2] and Chapter 4 of pj), is originated 
from the interaction of GCRs with the interplanetary magnetic field (IMF) and its irregularities. 
The IMF is the magnetic field that is carried outwards during the solar wind expansion. The 
interplanetary conditions vary as a function of the solar cycle which approximately lasts eleven 
years. In a solar cycle, when the maximum activity occurs, the IMF reverse his polarity. Thus, 
similar solar polarity conditions are found almost every 22 years [1]. In the HelMod Monte 
Carlo code version 1.5 (e.g., see Ref. [2]), the "classical" description of IMF, as proposed by 
Parker [S], was implemented together with the polar corrections of the solar magnetic field 
suggested subsequently in [51 Ej . This IMF was used inside the HelMod [2] code to investigate 
the solar modulation observed at Earth and to partially account for GCR latitudinal gradients, 
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i.e., those observed with the Ulysses spacecraft [HI [9]. In order to fully account for both the 
latitudinal gradients and latitudinal position of the proton-intensity minimum observed during 
the Ulysses fast scan in 1995, the HelMod Code was updated to the version 2.0 to include a 
new treatment of the parallel and diffusion coefficients following that one described in Ref. [lOj. 
In the present formulation, the parallel component of the diffusion tensor depends only on the 
radial distance from the Sun, while it is independent of solar latitude. 



2. The Interplanetary Magnetic Field 

Nowadays, we know that there is a Solar Wind plasma (SW) that permeates the interplanetary 
space and constitutes the interplanetary medium. In IMF models the magnetic-field lines are 
supposed to be embedded in the non-relativistic streaming particles of the SW, which carries 
the field with them into interplanetary space, producing the large scale structure of the IMF and 
the heliosphere. The "classical" description of the IMF was proposed originally by Parker (e.g., 
see [21 [5l [HI [121 [131 HI] and Chapter 4 of |H1). He assumed i) a constant solar rotation with 
angular velocity (w), ii) a simple spherically symmetric emission of the SW and iii) a constant 
(or approaching an almost constant) SW speed (T4w) at larger radial distances (r), e.g., for 
r > rb ^ IOR0 (where Rq is the Solar radius), since beyond rj, the wind speed varies slowly 
with the distance. The "classical" IMF can be analytically expressed as jl5] 



Bpar = ^ (er- - Te^) [1 -2H{e- 9')] , (1) 

where ^ is a coefficient that determines the IMF polarity and allows |Bpar| to be equal to -B9, 
i.e., the value of the IMF at Earth's orbit as extracted from NASA/GSFC's OMNI data set 
through OMNIWeb |16[ I17| : e^ and are unit vector components in the radial and azimuthal 
directions, respectively; 9 is the colatitude (polar angle); 9' is the polar angle determining 
the position of the Heliospheric Current Sheet (HCS)[18j; H is the Heaviside function: thus, 
[1 — 2H{9 — 9')] allows Bp^r to change sign in the two regions above and below the HCS [18j of 
the heliosphere; finally, 

^ ^ bjir — rh) sm9 , , 

r = tan ^' = — (2) 

''sw 

with ^ the spiral angle. In the present model u is assumed to be independent of the heliographic 
latitude and equal to the sidereal rotation at the Sun's equator. The magnitude of Parker field 
is thus: 

Bpar = ^\/l + r2. (3) 

In 1989 [6j, Jokipii and Kota have argued that the solar surface, where the feet of the field 
lines lie, is not a smooth surface, but a granular turbulent surface that keeps changing with 
time, especially in the polar regions. This turbulence may cause the footpoints of the polar field 
lines to wander randomly, creating transverse components in the field, thus causing temporal 





Period 


Years 






I) 


yl < Ascending 


1964.79- 


-1968.87, 


1986.70-1989.54, 2008.95-2009.95 


n) 


yl < Descending 


1964.53- 


-1964.79, 


1979.95-1986.70, 2000.28-2008.95 


III) 


A > Ascending 


1976.20- 


-1979.95, 


1996.37-2000.28 


IV) 


A > Descending 


1968.87- 


-1976.20, 


1989.54-1996.37 



Table 1. Definition of Ascending and Descending periods 
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Figure 1. (a) Maximum percentage difference between Bp^i and Bpar as a function of the 
solar distance inside the colatitude regions 20° < 6^ < 160° (620°), 25° < 6 < 155° (6250) and 
30° < 6 < 150° (6*30°). (b) Maximum value allowed of 5m as a function of colatitude: values 
in the "allowed" region guarantee that stream lines originated from the polar magnetic field do 
not cross the equatorial plane. Since Eq. ([5]) is symmetric with respect to the solar equatorial 
plane, values greater than 90° of colatitude lead to same results of those presented. 



deviations from the smooth Parker geometry. The net effect of this is a highly irregular and 
compressed field line. In other words, the magnitude of the mean magnetic field at the poles is 
greater than in the case of the smooth magnetic field of a pure Parker spiral. Jokipii and Kota [6] 
have therefore suggested that the Parker spiral field may be generalized by the introduction of a 
perturbation parameter [5{9)] which amplifies the field strength at large radial distances. With 
this modification the magnitude of IMF, Eq.(l3]), becomes [6]: 



Bpoi 



A 



i + r2 + 



5{ef 



(4) 



The difference of the IMF obtained from Eq. ^ and Eq. (jH) is less than ~1% for colatitudes 
20° < 6 < 160° (e.g, see Fig. [T^) and increases for colatitudes approaching the polar regions 
(e.g., see Figure 2 of [2]). 

In the present treatment, the heliosphere is divided into polar regions and a equatorial region 
where different description of IMF are applied. In the equatorial region the Parker's IMF, 
Eq. ([3]), is used, while in the polar regions we used a modified IMF that allows a magnitude as 
in Eq. (g]): 



B 
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B 



_ A 

Par — ~o 



r uj{r-rb) sine 
Br H d{9)eg — e, 



[1 -2H{e- 6')] Polar regions 



uj{r — Tb) sin( 



[1 -2H{9- e')] 



(5) 



Equatorial region. 



where equatorial regions are those with colatitude X° < 9 < (180° — X°). The symbol 9x° 
indicates the corresponding polar regions. 
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HelMod Energy bin KET channel 

(0.35 - 0.98) GeV (0.4 - 1.0) GeV 

(0.76 - 2.09) GeV (0.8 - 2.0) GeV 

(2.09 - 200) GeV > 2 GeV 

Table 2. Kinetic energy bins (in GeV) selected with HelMod Code and the corresponding 
proton energy channel for KET instruments on board the Ulysses spacecraft 



In order to have a divergence- free magnetic- field we require that the perturbation factor [d{9)] 
has to be: ^ 

^^^^ = [l-2g(g- noising ' 

where 5m is the minimum perturbation factor of the field. The perturbation parameter is let 
to grow with decreasing of the colatitude. However, in their original work, Jokipii and Kota [6] 
estimated the value of the parameter 6 between 10~^ and 3 x 10"'^. 

Since the polar field is only a perturbation of the Parker field, it is a reasonable assumption 
that stream lines of the magnetic field do not cross the equatorial plane, thus, remaining 
completely contained in the solar hemisphere of injection. This allows one to estimate an upper 
hmit on the possible values of 6m [see Fig. [TJb)]. Currently, we use 5m — 1 ^ 10 ^ by comparing 
simulations with observations at Earth orbit during Solar Cycle 23 (see Sect. [5]). 



Eneri 


gy range 


^lat,min 


(in dej 


;rees) 






^20° 


^25° 


6*30° 


0.35 - 
0.76 - 
2.09 - 


0.98 GeV 
2.09 GeV 
200 GeV 


-^tt - 
-5;? - 

- 


0+4 

-»_4 

_7+7 
' -7 


^-3 



Table 3. Latitudinal positions of mininiuin proton intensity (^iat,min) (i-'^ degrees) cis ci function 
of the kinetic energy, using three values for the extension of the polar regions. 



Eneri 


gy range 




An-s (%) 








^20° 


^25° 


^30° 


0.35 - 


0.98 GeV 




-lll^ 


-1311 


0.76 - 


2.09 GeV 






-9^1 


2.09 - 


200 GeV 




0+3 
"^-3 


-41^ 



Table 4. Percentages of North-South asymmetry of proton intensities (An-s) as a function of 
the kinetic energy, using three values for the extension of the polar regions. 



3. The Propagation Model 

Parker in 1965 [1] (see also Ref. [2], Chapter 4 of |3j and references therein) treated the 
propagation of GCRs trough the interplanetary space. He accounted for the so-called adiabatic 
energy losses, outward convection due to the SW and drift effects. In the heliocentric system 
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Enerj 


5y range 




Amax (%) 








^20° 


^25° 


^30° 


0.35 - 


0.98 GeV 


-34+5 


-35^1 


-36+3 


0.76 - 


2.09 GeV 


-22+^ 


-23t| 


-25+3 


2.09 - 


200 GeV 


-8i^ 


_q+3 





Table 5. Differences in percentage between the maximum and minimum proton intensities 
(Amax) as a function of the kinetic energy, using three values for the extension of the polar 
regions. 



the Parker equation is then expressed (e.g. see [2} [T9]): 



at - 8.. K^^g^ ) - ite-K''- + ""-'"l + slSf <'"='™' 

where U is the number density of particles per unit of particle kinetic energy T, at the time 
t- Vsw,i is the solar wind velocity along the axis Xj [20], Kfj is the symmetric part of diffusion 
tensor [1], Vd is the drift velocity that takes into account the drift of the particles due to the 
large scale structure of the magnetic field [211 ESI [23] and, finally, 

_ T + 2mrC^ 
I + m^c^ 

where rur is the rest mass of the GCR particle. The last term of Eq. ([7|) accounts for adiabatic 
energy losses [Il[23]. The number density U is related to the differential intensity J as ([2| I25j. 
Chapter 4 of [3] and references therein): 

where v is the speed of the GCR particle. 

Equation ([7|) was solved using the HelMod code (see the discussion in Ref. [2]). This treatment 
(i) follows that introduced in Refs. [Ml [23 [Ml [221 [33 [E] and (ii) determines the differential 
intensity of GCRs using a set of approximated stochastic differential equations (SDEs) which 
provides a solution equivalent to that from Eq. ([7|). The equivalence between the Parker equation, 
that is a Fokker-Planck type equation, and the SDEs is demonstrated in |28t I31j. In the present 
work, we use a 2D (radius and colatitude) approximation for the particle transport. The model 
includes the effects of solar activity during the propagation from the effective boundary of the 
heliosphere down to Earth's position. 
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Figure 2. Latitudinal relative intensity at r = lAU, obtained at different solar colatitudes for 
protons in the energy range defined in the Table [2] and using three definitions of polar regions: 
020°, 025° and O3Q0. 



The set of SDEs for the 2D approximation of Eq. ([7]) in heliocentric spherical coordinates is 
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Figure 3. For the polar regions 9 < 30° and 6 > 150°, latitudinal relative intensity at r =1 AU, 
accounting for GCR particles with 0° < 6' < 180° (^^o), 1° < f < 179° (6'r,i) and 2° < 6* < 178° 
(^R,2)) respectively, are shown as a function of the proton kinetic energy and solar colatitude. 
Results with 10° < 6 < 170° (6r,io) are comparable with those obtained with ^r,_2- 
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where fj,{9) = cos 6 and uji is a random number following a Gaussian distribution with a mean 
of zero and a standard deviation of one. The procedure for determining the SDEs can be found 

mm- 

In a coordinate system with one axis parallel to the average magnetic field and the other two 
perpendicular to this the symmetric part of the diffusion tensor Kf, is (see e.g. [32]): 




























(10) 



with A'li the diffusion coefficient describing the diffusion parallel to the average magnetic field, 
and K±^r and K±^0 are the diffusion coefficients describing the diffusion perpendicular to the 
average magnetic field in the radial and polar directions, respectively. In this work is that 
one proposed by Strauss and collaborators in Ref. [lOj (see also [331 1311 [35]): 

K\\ = ^Ko-^ (l + -4-) , (11) 



3 IGV V 1 AU 

where Kq is the diffusion parameter - described in Section 2.1 of [2j -, which depends on solar 
activity and polarity, /3 is the particle speed in unit of speed of light, P = pc/\Z\e is the particle 
rigidity expressed in GV and, finally, r is the heliocentric distance from the Sun in AU. 

In the current treatment, X|| has a radial dependence proportional to r, but no latitudinal 
dependence. Mc Donald and collaborators (see Ref. [36j) remarked that i) a spatial dependence 
of X|| - like the one proposed here - can affect the latitudinal gradients at high latitude and 
ii) it is consistent with that originally suggested in Ref. [6]. Furthermore, the perpendicular 
diffusion coefficient is taken to be proportional to K^^ with a ratio K}_^i/K\\ = 0.13 for both r 
and 6 i-coordinates. The latter value is discussed in Sect.0 

In addition, the practical relationship between Kq and monthly Smoothed Sunspot Numbers 
(SSN) [37j values - discussed in Section 2.1 of [2] - is currently updated using the most recent 
data from Ref. [35]. As in Ref. |2j, the Kq data are subdivided into four sets, i.e., ascending and 
descending phases for both negative and positive solar magnetic- field polarities (Table [T]). It has 
to be remarked that after each maximum the sign of the magnetic field (i.e., the A parameter 
in Eq. ([5])) is reversed. The updated practical relationships between Kq in AU^GV"^ s~^ and 
SSN values for 1.4 <SSN< 165 for the four periods (from I up to IV, listed in Table [1]) are 

I) Kq = 0.000297-2.9-10-6sSN+8.M0-9SSN2 + 1.46-10-10SSN3-8.4-10-13ssN'^ (12a) 

. / 0.000304846 - 5.8 • lO-^SSN if SSN <= 20 , . 

ii) ^0- I 000195 _ 2.3. 10-10SSN2 + 9. 1-10-5 if SSN > 20 ^^"^^^ 

III) Kq = 0.0002391 - 8.453 • lO-'^SSN (12c) 

IV) Kq = 0.000247 - 1.175 • lO^^SSN. (12d) 

The rms (root mean square) values of the percentage difference between values obtained with 
Eqs. (fl^ from those determined using the procedure discussed in Section 2.1 of [2] applied to the 
data from Ref. [38] were found to be 6.0%, 10.1%, 7.0% and 13.2% for the period I (ascending 
phase with ^ < 0), II (descending phase with A < 0), III (ascending phase with ^ > 0) and IV 
(descending phase with A > 0), respectively. 



4. The Magnetic Field in the Polar Regions 

Section [2] describes an IMF following the Parker Field with a small region around the poles 
in which such a field is modified. As already mentioned (see e.g.[6l |39l SQl HH US] US])) the 
correction is needed to better reproduce the complexity of the magnetic field in those regions. 
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5^ (xlO-^) 


Pk 


"R" model 


"L" model 


No Drift 


0.0 


0.10 


14.1 


11.0 


33.4 


1.0 


0.10 


11.7 


8.7 


33.2 


2.0 


0.10 


11.6 


8.3 


33.7 


3.0 


0.10 


11.6 


8.3 


33.7 


1.0 


0.11 


6.4 


9.0 


27.7 


2.0 


0.11 


7.8 


9.0 


28.3 


3.0 


0.11 


7.5 


8.8 


29.2 


1.0 


0.12 


6.3 


7.1 


23.5 


2.0 


0.12 


6.3 


7.3 


24.7 


3.0 


0.12 


7.1 


6.9 


24.4 


1.0 


0.13 


6.3 


6.4 


20.1 


2.0 


0.13 


6.6 


7.6 


20.4 


3.0 


0.13 


6.7 


7.7 


20.5 


1.0 


0.14 


7.3 


7.0 


15.9 


2.0 


0.14 


7.3 


7.2 


16.4 


3.0 


0.14 


7.2 


6.5 


16.8 



Table 6. Average values (last three columns) of ?7rins (in percentage, %) as a function of 
5m (xlO-5) and pk, for BESS-1997, AMS-1998, PAMELA-2006/08, obtained from Eq. ([16]) 
without enhancement of the diffusion tensor along the polar direction [pE = 1), using "R" and 
"L" models for the tilt angle and No Drift approximation. The differential intensities were 
calculated accounting for particles inside the heliospheric regions for which solar latitudes are 
lower than 1 5.7° I . 



Moreover, Ulysses spacecraft (see e.g. [HI SSI SB] ) explored the heliosphere outside the ecliptic 
plane up to ±80° of solar latitude at a solar distance from ~ 1 up to ~ 5AU. Using these 
observations, the presence of latitudinal gradient in the proton intensity could be determined 
(e.g., see Figure 2 of Ref. |9j and Figure 5 of Ref. [8j). The data collected during the latitudinal 
fast scan (from September 1994 up to August 1995) show (a) a nearly symmetric latitudinal 
gradient with the minimum near ecliptic plane, (b) a southward shift of the minimum and 
(c) an intensity in the North polar region at 80° exceeding the South polar intensity. In Ref. 
[9] a latitudinal gradient of ~ 0.3%/degree for proton with kinetic energy > 0.1 GeV was 
estimated. While in Ref. [8j the analysis to higher energy was extended estimating a gradient 
of ~ 0.22%/degree for proton with kinetic energy > 2 GeV. The minimum in the charged 
particle intensity separating the two hemispheres of the heliosphere occurs ~ 10° South of 
the heliographic equator [9]. In addition, an independent analysis that takes into account the 
latitudinal motion of the Earth and IMPS confirms a significant (~ 8° it 2°) southward offset of 
the intensity minimum [9j for T > 100 MeV proton. Furthermore in Ref. [8], a southward offset 
of about ~ 7° is evaluated; this offset of the intensity minimum results to be independent of the 
particle energy up to 2 GeV. Finally, in Ref. [9], the intensity in the North polar region at 80° 
is observed to exceed the South polar intensity of ~ 6% for protons with T > 100 MeV. 

Using the present HelMod code (version 2.0), we could investigate i) the latitudinal gradient 
of GCR intensities resulting from solar modulation and ii) how the magnetic-field structure of 
the polar regions, as defined in Sect. [21 is able to influence the GCR spectra on the ecliptic plane. 
As previously defined, we denote with 9x° a polar region of amplitude X° from polar axis, i.e., 
6* < X° and 6* > 180° - X°. Three regions with X° = 20°, 25° and 30° were investigated. 
Outside any of these regions, the ratio between Bpoi and Bpar in Eq. ([5]) is less than ~1% (see 
Fig. [TJ and, thus, it ensures a smooth transition between polar and equatorial regions. For the 
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^ I I I I I I I I I I I I 

10"' 10° 10' [GeV] 

Kinetic energy 

Figure 4. Proton differential intensity determined with the HelMod code (continuous line) 
compared to the experimental data of BESS-1997; the dashed line is the LIS (see the text). 



purpose or this study we consider an energy binning closer to those presented in Ref. [8]. The 
KET instrument [U] collects proton data in three "channels" one with energies ranging from 
0.038 GeV up to 2.0 GeV and two for particle with kinetic energy T > 0.1 GeV and T > 2 GeV, 
respectively. A successive re-analysis of the collected data allowed the authors to subdivide 
the 0.25-2 GeV "channel" in three "sub-channels" of intermediate energies. Since the Present 
Model is optimized - as discussed in Ref. [2] - for particles with rigidity greater than 1 GV 
(i.e., ~ 0.444 GeV), the present results are compared only with the corresponding "channel" or 
"sub-channels" suited for the corresponding energy range (see Table [2]). 

At 1 AU and as a function of the solar colatitude, the GCR intensities for protons are shown 
in Fig. [21 For a comparison with Ulysses observations, the modulated intensities of protons - 
resulting from HelMod code - were investigated from 80° (North) and down to —80° (South). 
They were obtained using the HelMod code and selected using the energy bins reported in table 
[21 In Fig. [21 the latitudinal intensity distribution is normalized to the corresponding South Pole 
intensity. The quoted errors include statistical and systematic errors. The distributions were 
interpolated using a parabolic function expressed as: 

Ii9i.t) = a + ci9i,t + df, (13) 

where /(^lat) is the normalized intensity, ^lat is the latitudinal angl^l] and a,b and c are 
parameters determined from the fitting procedure. The so obtained fitted curves are shown as 
continuous lines in Fig. [21 Furthermore, the latitudinal positions of niininiuni intensity (^lat min); 



^ The latitudinal angle is 6'iat = 90° - 0. 
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Figure 5. Proton differential intensity determined with the HelMod code (continuous hne) 
compared to the experimental data of AMS-1998; the dashed line is the LIS (see the text). 



percentages of North-South asymmetry of intensities (An-s) and differences in percentage 
between the maximum and minimum intensities (Amax) were also determined from the fitting 
procedure and are listed in Table O Table S] and Table [5l respectively. The quoted errors - 
following a procedure discussed in Ref. - are derived varying the fitted parameters in order 
to obtain a value of the parameter ryrms,iat two times larger than the one resulting from the best 
fit (/best)- ??rms,iat IS defined as: 



with 



^rms.lat 



— /best(^j) 



(14) 



(15) 



est \ 



where 6i is the central value of the ith. latitudinal bin of the differential intensity distribution 
and crn,i,iat are the errors including the experimental and Monte Carlo uncertainties. For ^30°, 
i.e., assuming a modified polar magnetic-field for 6 < 30° and 9 > 150°, we found a general 
agreement with Ulysses observations. The position of ^lat.min is compatible within the errors 
with one observed in Ref. [8j, as well as the values of An-s and Amax- 

The HelMod code allows one to investigate the relevance of the treatment of the polar region 
magnetic-field with respect to the resulting modulation of GCR. Thus, the latitudinal normalized 
intensities were obtained excluding a few (small) regions nearby the poles. This was determined 
from reducing the latitudinal spatial phase-space admissible for pseudo-particles (see Ref. [2]) 
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Figure 6. Proton differential intensity determined with the HelMod code (continuous hne) 
compared to the experimental data of PAMELA 2006/08; the dashed line is the LIS (see the 
text). 



- i.e., the latitudinal extension of GCR particles taken into account - to 1° < 9 < 179° (^r,i), 
2° < 6* < 178° (6Ir,2) and 10° < 6 < 170° (6Ir,io)- The so obtained latitudinal gradients are 
compared with the full latitudinal extension, (0° < 6* < 180°), in Fig. [3l By an inspection 
of Fig. [3l one may lead to the conclusion that the GCR diffusion nearby the polar axis has a 
large impact on the latitudinal gradients in the inner heliosphere. As a consequence, the IMF 
description in the polar regions is relevant in order to reproduce the observed modulated GCR 
spectra. 

5. Comparison with Observations During Solar Cycle 23 

The agreement of HelMod simulated spectra with observations during solar cycle 23 is 
investigated via quantitative comparisons using Eqs. and (fT7|) . However, since the structure 
of the heliosphere is different in high and low solar activity the two periods are separately 
analyzed. 

The HelMod Code ^ (version 2.0) allowed us to investigate how the modulated (simulated) 
differential intensities are affected by the (1) particle drift effect, (2) polar enhancement of the 
diffusion tensor along the polar direction {K^^q, e.g., see Ref. [IS]), and, finally, (3) values of the 
tilt angle (at) calculated following the approach of the "R" and "L" models [39]. This analysis 
also allow us to estimate the values of IMF parameters that better describe the modulation 
along the entire solar cycle. The effects related to particle drift were investigated via the 
suppression of the drift velocity (No Drift approximation), this accounts for the hypothesis 
that magnetic drift convection is almost completely suppressed during solar maxima. The 
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5^ (xlO- 


Pk 


"R" model 


"L" model 


No Drift 


0.0 


0.10 


11.2 
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15.4 
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11.0 


10.1 


15.8 


2.0 


0.10 


9.6 


10.0 


16.7 


3.0 


0.10 


9.6 


10.0 


16.7 


1.0 


0.11 


13.4 


13.1 


16.0 


2.0 


0.11 


12.7 


12.9 


15.4 


3.0 


0.11 


12.7 


12.5 


16.2 


1.0 


0.12 


1 O "7 

18.7 


1 '1 '1 

17.7 


13.4 


2.0 


0.12 


18.3 


16.9 


12.8 


3.0 


0.12 


18.1 


17.3 


12.8 


1.0 


0.13 


23.3 


23.5 


14.3 


2.0 


0.13 


25.0 


24.7 


13.3 


3.0 


0.13 


24.3 


24.2 


13.1 


1.0 


0.14 


32.3 


30.7 


18.0 


2.0 


0.14 


32.8 


30.8 


17.1 


3.0 


0.14 


31.5 


30.7 


17.9 



Table 7. Average values (last three columns) of ?7rins (in percentage, %) as a function of 
5m (xlO-^) and pk, for BESS-1999, BESS~2000, BESS-2002, obtained from Eq. ^ without 
enhancement of the diffusion tensor along the polar direction [pE = 1), using "R" and "L" 
models for the tilt angle and No Drift approximation. The differential intensities were calculated 
accounting for particles inside the heliospheric regions for which solar latitudes are lower than 
|5.7°|. 



differential intensities were calculated for K^^^ = pEK±^r with values of pE of 1, 8 and 10, 
i.e., no enhancement, that suggested in Ref. [50] and that suggested in Ref. [IHl [2] (and 
reference there in), respectively. Furthermore, the modulated proton spectra were derived from 
a LIS whose normalization constant depends on the experimental set of data and were already 
discussed in Ref. [2]). In addition, the differential intensities were calculated accounting for 
particles inside heliospheric regions where solar latitudes are lower than |5.7°|. 

During the period of high solar activity for the solar cycle 23, the BESS collaboration took 
data in the years 1999, 2000, and 2002 (see sets of data in Ref. [51]). For period not dominated 
by high solar activity in solar cycle 23, BESS, AMS and PAMELA collaborations took data, 
i.e., BESS-1997 [ST], AMS-1998 P], and PAMELA-2006/08 [53]. 

Following the procedure described in Ref. [2j, the observation data were compared with those 
obtained from HelMod code using the error-weighted root mean square (r/rms) of the relative 
difference (??) between experimental data (/cxp) and those resulting from simulated differential 
intensities (/sim)- For each set of experimental data and with the approximations and/or models 
described above, we determined the quantity 

with 

fsim{Ti) — /exp(2~i) 

Vi = r~rF\ ' (17) 

where Tj is the average energy of the ith energy bin of the differential intensity distribution and 
arj^i are the errors including the experimental and Monte Carlo uncertainties; the latter account 
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Figure 7. Proton differential intensity determined with the HelMod code (continuous hne) 
compared to the experimental data of BESS-1999; the dashed line is the LIS (see the text). 



for the Poisson error of each energy bin. The simulated differential intensities are interpolated 
with a cubic spline function. The modulation results are studied varying the parameters 6'^^ - 
from 0, i.e., the non-modified Parker IMF, up to 3 x 10~^, see Fig. mb) -, K±^r/I^\\ = Pk (from 
0.10 up to 0.14) and K^^/K^^j. = pE (from 1 up to 10) seeking a set of parameters set that 
minimize r/rms- In Tabled the average values of r/rms (in percentage, %) for low solar activity 
periods are listed. They were obtained in the energy range □ from 444 MeV up to 30 GeV using 
the "L" and "R" models for the tilt angle at and for the No Drift approximation and without 
any enhancement of the diffusion tensor along the polar direction {K±^^). The results derived 
with the enhancement of the diffusion tensor along the polar direction indicate that for pe = 8 
and 10 one obtains a value of r/rms that is from 1.5 up to 3 times larger with respect the case 
without enhancemenllfl. From inspection of Table El one can remark that the drift mechanism 
leads to a better agreement with experimental data. Furthermore the, "R" and "L" models 
for tilt angles are comparable within the precision of the method (discussed in Ref. [2]). The 
minimum difference with respect to the experimental data occurs when = 0.11 — 0.13 and 
6m = 1.0 X 10"^ for both "R" and "L" models, with the "L" model shghtly preferred to "R". 

In Figs. HHSl the differential intensities determined with the HelMod code are shown 
and compared to the experimental data of BESS-1997, AMS-1998, and PAMELA-2006/08, 
respectively; in these figures, the dashed lines are the LIS as discussed in Ref. |2j. These 
modulated intensities are the ones calculated for a heliospheric region where solar latitudes are 

^ Above 30 GeV, the differential intensity is marginally (if at all) affected by modulation. 

^ For a comparison, the scalar approximation presented in Ref. [2], i.e., assuming that the diffusion propagation 
is independent of magnetic structure, leads to and average T^rms of ~ 15%. 
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Figure 8. Proton differential intensity determined with the HelMod code (continuous hne) 
compared to the experimental data of BESS-2000; the dashed line is the LIS (see the text). 



lower than |5.7°|, using /o^ = 0.13, 5m = 1-0 x 10"^ and pE = 1 with the "L" model. Finally, one 
can remark that the present code, combining diffusion and drift mechanisms, is also suited to 
describe the modulation effect in periods when the solar activity is no longer at the maximum. 

In Table [7] we present the averages r/rms (in percentage, %) during the periods dominated by 
high solar activity. The simulated differential intensities were obtained for a heliospheric region 
where solar latitudes are lower than |5.7°| without any enhancement of the diffusion tensor along 
the polar direction (i^j,,^). The simulations with pE = 8 and 10 lead to 77rnis comparable with 
those presented in Table [71 However, using pE = 1 provides a better agreement to experimental 
data at lower energy. From inspection of Table \7\ one can note that "R" and "L" models for tilt 
angles yield comparable results within the precision of the method. Furthermore the minimum 
difference with the experimental data occurs when pk = 0.10 and 6^ = (2.0-3.0) x 10"^ with 
the"L" model slightly preferred to "R". The HelMod parameter configuration, which minimizes 
the difference to the experimental data, are reported in Table [H] one may remark that the 
No Drift approximation is (almost) comparable to a drift treatment for both BESS-2000 and 
BESS-2002 with data collected during and after the maximum of the solar activity. Apparently, 
the drift treatment is needed in order to describe BESS-1999 with data taken during a period 
approaching the solar maximum. 

In Figs. ©-dHl), the differential intensities determined with the HelMod code are shown and 
compared with the experimental data of BESS-1999, BESS-2000, and BESS-2002, respectively; 
in these figures, the dashed lines are the LIS as discussed in Ref. [2]. These modulated intensities 
are the ones calculated for a heliospheric region where solar latitudes are lower than |5.7°|, using 
K±^fj^ = -ftr_L,r independently of the latitude and including particle drift effects with the values 



- 16 - 




Kinetic energy 



Figure 9. Proton differential intensity determined with the HelMod code (continuous hne) 
compared to the experimental data of BESS-2002; the dashed line is the LIS (see the text). 



Observations 


"R" model 


"L" model 


No Drift 


BESS-1999 


9.3 


10.6 


25.7 


BESS-2000 


12.5 


12.6 


16.7 


BESS-2002 


6.9 


6.7 


7.7 



Table 8. Average r/^ms (in percentage, %), for BESS-1999, BESS-2000, BESS-2002, obtained 
from Eq. ()16p without enhancement of the diffusion tensor along the polar direction (pE = 1), 
6m = 2.0 X 10-^ pk = 0.10 and using "R" and "L" models for the tilt angle and No Drift 
approximation. The differential intensities were calculated accounting for particles inside the 
heliospheric regions for which solar latitudes are lower than |5.7°|. 



of the tilt angle from the "L" model. Finally, it is concluded that the present code combining 
diffusion and drift mechanisms is suited to describe the modulation effect in periods with high 
solar activity [21 [501 [Si]- 

6. Conclusion 

In this work an IMF, which combines the Parker Field and its polar modification, is presented. 
In the polar regions, the Parker IMF was modified with an additional latitudinal components 
according to those proposed by Jokipii and Kota in Ref. [6j. We found the maximum perturbed 
value with this component yielding, as a physical result, streaming lines completely confined in 
the solar hemisphere of injection. 

The proposed IMF is, then, used within the HelMod Monte Carlo code to determine the 
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efFects on the differential intensity of protons at lAU as a function of the extension of polar 
region, in which the modified magnetic-field is employed. We found that a polar region contained 
within 30° of colatitudc is that one ensuring a very smooth transition to the equatorial region and 
allows to reproduce qualitatively and quantitatively the latitudinal profile of the GCR intensity, 
and the latitudinal dip shift with respect to the ecliptic plane. Finally we determined how the 
polar region diffusion is mostly responsible of the proton intensity latitudinal gradient observed 
in the inner heliosphere with the Ulysses spacecraft during 1995. 
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